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Abstract. We show that choreographic three bodies {x(t),x(i -I- T/3),x{t — T/3)} 
of period T on the lemniscate, x(t) = (x -I- ycn(t))sn(i)/(l -I- cn^(t)) parameterized by 
the Jacobi's elliptic functions sn and cn with modulus k'^ — {2 + \/3)/4, conserve the 
center of mass and the angular momentum, where x and y are the orthogonal unit 
vectors defining the plane of the motion. They also conserve the moment of inertia, the 
kinetic energy, the sum of square of the curvature, the product of distance and the sum 
of square of distance between bodies. We find that they satisfy the equation of motion 
under the potential energy I]i<j (1/2 - \/3/24rfj) or J2i<j 1/2 ^^nj - J2i VS/Srf, 
where the distance between the body i and j, and the distance from the origin. 
The first term of the potential energies is the Newton's gravity in two dimensions 
but the second term is the mutual repulsive force or a repulsive force from the 
origin, respectively. Then, geometric construction methods for the positions of the 
choreographic three bodies are given. 



PACS numbers: 45.20.Dd, 45.50.Jf 



1. Introduction 



Choreographic motion of bodies is a periodic motion on a closed orbit, bodies 
chase each other on this orbit with equal time-spacing. Recently, choreographic motions 
under the Newton's gravity are found and paid attentions. C. Moore found a figure 
eight three-body choreographic solution by numerical calculations. A. Chenciner and 
R. Montgomery gave a rigorous proof of the existence of the choreographic figure 
eight three-body solution. At the same time, C. Sim6|^, ^ found many remarkable 
choreographic A^-body solutions by numerical calculations. Although exact form of these 
choreographic A^-body solutions are still unknown, Simo's figure eight choreographic 
three-body solution is very similar to an affine transformed lemniscate 

Concerning the relation between the lemniscate and the Newton's equation of 
motion, it is well known |P] that a point particle on the lemniscate = cos (26') with 



Choreographic Three Bodies on the Lemniscate 



2 



6 = 1/2 sin~^(2/t), —1 < 2lt < 1 satisfies the equation of motion under tlie central 
potential U{r) = —P/{2r^). Here r, 6, I, and t represent the radius, the azimuthal 
angle, the angular momentum of the particle and the time, respectively. This motion 
is not periodic. This particle starts from the origin at t = —1/(2/), travels right leaf 
of the lemniscate, and finally collides with the origin at t = 1/(2Z). Two-body problem 
on the lemniscate is derived from the one body problem. Two particles of equal masses 
start from the origin at t = —1/(2/) to opposite directions, each particle travels the 
left leaf or right leaf of the half size lemniscate, and collides with each other at the 
origin at t = 1/(2/). No analytic solution is known about more than three bodies on 
the lemniscate. 

Being stimulated by the Moore, Chenciner, Montgomery and Simo's remarkable 
study, we investigated physical and geometrical properties of the choreographic motion 
of three bodies on the lemniscate. We describe the results in section 2 and proof is 
given in section ^. In section 4, we discuss the relation between the choreographic three 
bodies on the lemniscate and the rectangular hyperbola, and give a geometrical method 
to determine the positions of the choreographic three bodies. Section 5 is the summary. 



2. Choreographic three bodies on the lemniscate 

Let us parameterize by using the Jacobi's elliptic functions sn and cn the lemniscate 
x(t) = x(t)x + y{t)y which satisfies (x^ + y^)^ = — y"^ as follows, 

sn(t) 

^^'^ = TwM ' 

_ sn(t)cn(t) ^ ^ 

' " l + cn2(t)' 

where x and y are the two basic orthogonal unit vectors defining the plane of the motion. 
This is a smooth periodic motion on the lemniscate with period 

T = AK, (2) 

where K is the complete elliptic integral of the first kind. Positions of the choreographic 
three bodies are 

{xi (t) , X2 (t) , X3 (t) } = {x(t) , x(t + 4K/3) , x(t - 4A73) } . (3) 

In the following, we use notations v(t) = dx(t)/dt and a(t) = d^x(t)/dt^. 
Straight calculation shows a relation between the curvature 
1 _ |v(^) X a(t)| 

|v(t)|3 

and the distance from the origin, 

p~'{t) = 9x2(t), (4) 

which is, of course, parameterization invariant relation. And also, we get a relation 
between the velocity and the distance from the origin, 

v^W + ik' - i)x2(t) = (5) 
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Figure 1. The lemniscate and the positions of x(i) with modulus fc^ = (2 + \/3)/4 at 
t = jK/3 for j = 0, 1,2, 11. Full circles • labeled by j represent the positions of 
x(i). 



for arbitrary modulus k. 

Conservation of the center of mass 

x(t) + x(t + 4if/3) + x(t - AK/3) = 

is satisfied if and only if the value of the modulus is 

J,2 2 + V/3 



(6) 



(7) 



as shown in section 0. Parameterization in (0), (^) and defines the motion of 
the choreographic three bodies on the lemniscate. In figure |I], the lemniscate and the 
positions x(t) at t = jK/3 for j = 0, 1, 2, 11 are shown. The order of magnitude of 
the label j represents the direction of the motion and the points having the same label 
in modulo 4 are the positions of the choreographic three bodies on the lemniscate 
at t = 0,K/3,2K/3,K,.... 

We find that this motion conserves the moment of inertia and angular momentum 



Ex? 

^ Xj X V, 



V3, 



0, 



as shown in section]^. Conservation of the moment of inertia and relations 
the conservation of the sum of square of the curvature and the kinetic energy, 

Ep7' = QV3, 



Simple algebra shows that 
conservation, 



(8) 
(9) 

(I) yield 

(10) 
(11) 



[) implies J2 
3V3. 



3J2i^'i- Thus we get another 

(12) 



i<j 
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The affine transformed motion x(t)x + k'^y{t)y is numerically close to the Simo's 
figure eight motion. It is difficult to distinguish the two orbits if they are plotted on the 
usual computer display. This similarity and the above conservation laws suggest that the 
choreographic three bodies on the lemniscate (|1]), (|^) and (|^ may satisfy an equation 
of motion under a interaction potential energy similar to that of the Newton's gravity. 
We find that it surely satisfies an equation of motion under the Newton's gravity in two 
dimensions with extra repulsive force, 
d^ 

_ 1 [ x(t + 4ir/3) - x(t) x(t - 4ir/3) - x(t) ] (^3) 

Newton ^ | ^^^^ ^ ^^^^^ _ + ^^^^ _ ^^^^^ _ j , 

as shown in section |. The repulsive force Fj-epuisive can be expressed in two ways 

Frcpulsivel = — (14) 

^/U + ^?)-xwUfx«-^^)-x«m, (15) 



■T repulsive2 ^2 1 V ' 3 ' --\-/ J ' g 

Due to the conservation of the center of mass, the two repulsive force is equal, 
Frepuisivei = Fi.epuisive2 5 ou the orbit. The corresponding potential energy for the equation 

of motion with Frepuisivei IS 

t^ = E^nr.,-E^x^ (16) 
and with Frepuisive2 is 

2 ^ 24 



y = i:\>r,-^rl\, (17) 



^2 



^3) 



where is the distance between the point i and j, i.e., = y (xj — x 

With the conservation of the kinetic energy and the moment of inertia, the 
conservation of total energy with potential energy f/, in (|16|) , implies the conservation 
of the product of rfy Evaluating the value at t = 0, we get 

rUtyutyut) = (18) 

3. Proof 

3.1. Conservation of the center of mass 

Since ^{K) = x, the x-component of the other two points must be x{K±4K/3) = —1/2. 
This equation gives, 

sn{K/3) = VS-l. (19) 

Evaluating sn(10-ft'/3) in two ways, 

sn(10K/3) = sn{3K + K/3) = -cn{K/3)/dn{K/3) 
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Table 1. Values of sn(t), cn(i), dn(i). 



t 


sn(t) 




dn(t) 


K/S 


\/3- 1 


3i/4(\/3- l)/V2 


1/V2 


2K/3 


3i/4(V3_ 1) 




(V3-l)/2 


4:K/3 


3i/4(V3_ 1) 


-2 + V3 


(\/3-l)/2 


5K/3 




-3i/4(\/3- 1)/V2 


1/V2 



and 



we get 



sn(10K/3) = sn(4ir - 2K/3) = \/ Jp^^/^^/gl) ^ 



2 1 - 2sn(ir/3) 



^ sn4(ir/3)-2sn3(ir/3)' 

Substituting the value of (|T9|) into (pOD gives the value of fc^ in (|^). Inversely, the value 
of the modulus in (|^ gives the value of sn(_ft'/3) in (|T9D . The value of sn(t), cn(t), dn(t) 
at some t are shown in table |l]. To complete the proof of (|^), it is convenient to use the 
complex variable 

and consider the complex t-plane. We take the fundamental cell for t-plane with 

-2K <m< 2K, -2K' <^t< 2K\ (21) 

where K' is the complementary elliptic integral of the first kind. Since x''^\t) has 4 
simple zeros in the fundamental cell at t = —2K — 2ii^', — 2ii^', —2K and 0, degree of 
the elliptic functions x^^\t) is 4. Thus they should have 4 poles in the fundamental cell. 
Note that, 

1 



X 



ksn(u) — dn(M) 

From the value of fc^ in (|^) and the table 0, at m = K/3 and 5K/3, ksn{u) = dn(M). 
So, x'-+\t) has poles at t = K/3 + iK' and hK/3 + iK' . Since = it 

also has poles at t = —{K/3 + iK') and —{bK/3 + iK'). We find 4 poles for degree 4 
elliptic function then all these poles are simple poles. Table || shows the poles 

and residues for x^^\t) = sn(t)/(l — i cn(t)). Here we used the following notation, 

ai = -K + iK', a2 = K/3 + iK' , = 5K/3 + iK'. (22) 

It is clear from the table ^that all these poles are cancelled in x^~^\t) + x^^\t + AK / 3) + 
x^^\t — 4:K/3). Therefore, it is constant. Evaluating this constant at t = 0, we get 

x^+^ (t) + (t + AK/3) + x^+^ {t - AK/3) = 0, (23) 

which is equivalent to (0). 
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Table 2. Poles and residues of sn(f)/(l — i cn(^)) and 1/(1 — i cn{t)). 



function 


t = a2 


t = as 


t = —a2 


t = -a3 


sn(i)/(l - i cn(i)) 


72/31/4 


-72/31/4 


v/2/31/4 


-V2/31/4 


l/(l-icn(i)) 


1/31/4 


-1/31/4 


-1/31/4 


1/31/4 



3.2. Conservation of the moment of inertia and the angular momentum 



Let us consider a function 
1 



Since 1/(1 — i cn(M -|- iK')) = ksn{u) / {ksn{u) — dn(-u)), the function (p^ has simple 
poles at t = ±02, ±^3 in the above cell (0). In table we list the poles and residues 
of this function. From this table and the value at t = 0, it is clear that the function 
(15) satisfies 

1 . ' + , ^ ' (25) 



1 - i cn(t) 1 - i cn(t + AK/i) 1 - i cn(t - 4K/3) 2 
Now, consider a function j^^\t) = x^^\t)dx^^\t) / dt. Since x*^^^ = x ±iy, 
■ ca-w X , dx dy. ., dy dx, 

On the other hand, using x'^^^(t) = sn(t)/(l =F i cn(t)), we get 

dt 1 — 1 cn(t) 



then, from equation (PSQ 

^■(+)(t) + ji+)(^t + 4K/3) + /+\t - 4K/3) = 0. (28) 
Equations ( p6|) and ( ^8]) prove the conservation of the moment of inertia and equation 



(I). We get equation (||) from the value at t = 0. 

3.3. Equation of motion 
Let us consider a function 

Ax^'\t) = x^-\t + AK/3) - x^-\t). 

The function x^^\t) has 4 simple poles at iag) ic^S) ^he star being complex conjugation. 
Then, the function Ax*-^-*(t) has 6 simple poles at ±a]|, ±027 ic^a- So, the degree of this 
function is 6. 

On t = u + iK', the function Ax^^\t) is 

Ax^^Hu + iK') = - - 

^ ^ A;sn(M + AK/3) + dn{u + AK/3) ksn{u) + dn(n) ' 

And from the table 0, we see sn(5-ft'/3) = sn(K/3) and dn(5K/3) = dn(i^/3). Thus, 
at M = K/3, i.e., t = K/3 + iK', Ax^~\t) has zero point. To see the behaviour of the 
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function around the zero point, let take u = K/3 + Au. Then, using sn{2K — z) = sn{z) 
and dia{2K — z) = dn(2;), we get 

/\x^~\K/?, + Au + iK') = -Ax^-\K/3 - Au + iK'), 

i.e., Ax^^\t) is an odd function around the zero point t = a2 = K/3 + iK'. 

Series expansion of x^^\u + AK/S + iK') and x^^\u + iK') around u = K/3 yields 



X 



(-), 



u + 4K/3 + iK') 



X 



^-\u + iK') 



0I/4 q o3/4 

Thus, we get 

0I/4 o3/4 

Ax(-)(t) = ^(t - a,)' + - a,f + ■ ■ • , (31) 

and 

Ax(")(t) 3i/4(t_«2)3 V2(t-a2) ' ^ ^ 

i.e., l/Ax*-~-'(t) has triple pole at t = a2- Since x^~\—t) = —x^^\t), Ax^^\t) also has 
triple zero at t = —0:3 = —5K/3 — iK', as follows, 

Aa;(^)(-a3 + At) = Ax^~\a2 - At) 

Thus 1/Ax^~\t) has triple pole at t = —a^ with principal part 



Ax(-)(t) 31/4 (t + «3)3 ^2 (t + «3) 

We found two triple poles for the elliptic function of degree 6. 

In table |, we list the principal part of l/Ax'^'^t), -l/Ax'^-\t - AK/3), 
d^x''+-*(t)/dt^ and x^~^\t). From this table, we see the following equation is satisfied. 

This is the equation of motion in (|1^) with the repulsive force in (p^. 
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Table 3. Principal part around the pole with a = l^/lji'^l^, b = 'i^l^j-^. 

function t = t = «3 t = —a^ t = —as 

^^"^ ^'''i {t-a2f t-a2 (t+as)^ ^ i+Q3 

-\AT^-^(f - ^]-^ k-^ ^ 

yb 3 J (t-Q3)3 ^ t-a:i (t+a2)^ t+a2 

d^X^+\t)/de a a a a 



it-a2)^ (t-aa)^ (t+a2)3 it+a'i)^ 

(+)(t) I 1 I ^ 

^ ' b{t-a2) b{t-a:i) b{t+a2) 6(i+a3) 



4. Choreographic three bodies and the rectangular hyperbola 

In general three-body problem on the plane, not restricted to this lemniscate case, the 
conservation of the center of mass and zero angular momentum has simple geometrical 
meaning: three tangent lines from the three bodies must meet at a point (infinity is 
allowed) at each instant. That is, there exist scalars Xi{t) and a vector c(t), such that, 

c = Xj + AjVj, for i = 1,2, 3. (35) 

The solution Aj is, 

A,= «^pi^i^^,withV,^i, (36) 

(V, X Vj) ■ z 



and the vector c(t) is given by 

_ liVj - IjWi 



[yi X Vjj ■ z 



, with (^,j) = (1,2), (2, 3), (3,1), (37) 



where z = x x y and Zj = (xj x Vj) ■ z. 

For the lemniscate, we find that the orbit of c(t) is the rectangular hyperbola as 
shown in figure |], 

(c ■ x)2 - (c ■ y)2 = 1 (38) 

by a direct algebraic calculation. 

Let us introduce terms "forward" and "backward" for the direction of the 
choreographic motion on the lemniscate. We call the motion "forward" if points pass 
through the origin upward (from lower left to upper right or from lower right to upper 
left) and call "backward" otherwise. In other words, a "forward" motion is a clockwise 
motion on the right leaf of the lemniscate or a anti-clockwise motion on the left leaf. 

We observed the following two properties, (i) For the "forward" motion of the 
choreographic three bodies, the point c always moves upward on the rectangular 
hyperbola, (ii) Four positions of the choreographic three bodies and the cross point 
on the rectangular hyperbola, Xj, i = 1,2,3 and c, are in the different quadrants each 
other as shown in figure]^. Notice that when the point c jumps the leaf of the hyperbola 
one of the three bodies passes through the origin. And when the point c passes through 
the horizontal axis upward (downward) one of the three bodies passes through the 
horizontal axis at the same point in opposite direction, downward (upward). 
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1.5 




Figure 2. Snapshot at t = —K/6. Three tangent Unes from the choreographie three 
bodies meet at a point on the rectangular hyperbola. Full circles • are positions of 
the three bodies. 

Then, we will have two questions as to how the choreographic three points can 
be determined geometrically. The first question: How can we find the positions of the 
choreographic three bodies for an arbitrarily given point c on the rectangular hyperbola? 
We can draw four (in general) tangent lines to the lemniscate from the point c. Therefore 
we have four (in general) contact points instead of three. How can we select the 
choreographic three points from among them? We observed that the following two 
methods work, (i) Move the point c slightly upward. Then, we will observe that three 
of the four contact points move "forward" , and one "backward" . Of course the former 
is the positions of the choreographic three bodies, (ii) Choose three points which is not 
in the same quadrants the given point c is in. 

The second question: How can we find other two choreographic positions for an 
arbitrarily given point xi on the lemniscate? The tangent line of the lemniscate which 
contact at the given point xi have two cross point, di, d2, on the rectangular hyperbola. 
How can we select one? Again we have two method, (i) If we move the point Xi 
"forward", one of dj will move upward, and the other downward. The former is the 
point c the point xi corresponds to. (ii) Select a cross point dj which is in the different 
quadrants from the one the given point Xi is in. Using the method described in the 
above paragraph, we can find positions X2 and X3 for the other two choreographic bodies. 

The above procedures give a geometrical method to find a set of positions of the 
choreographic three bodies on the lemniscate. 
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5. Summary 

We have shown that the choreographic three bodies on the lemniscate (|I]), (|^) and (|^) 
satisfies the equation of motion under two different potential energies, U defined in (^) 
and V in (p!?]). It means that the choreographic three bodies on the lemniscate can 
be realized by the following two ways, (i) Consider four particles, the infinitely heavy 
particle and three particles 1-3 of equal masses. Take the Newton's gravity in two 
dimensions 

Uij = ^Inrij (39) 

as the interaction potential energy between particle i and j with i,j = 1,2,3, and the 
artificial repulsive potential energy 

Uo^ = -^x? (40) 

as the interaction potential energy between the particle and i with i = 1,2, 3. Then, 
put the infinitely heavy particle at the origin and put the three particles 1-3 on the 
lemniscate according to equation (^. (ii) Consider three particles 1-3 of equal masses. 
Take the Newton's gravity in two dimensions accompanied with the artificial repulsive 
potential 



V^, = l^nr,,-^rl (41) 



as the interaction potential energy between particle i and j with i,j = 1,2,3. Then, 
put the three particles 1-3 on the lemniscate according to equation (^. 

As a consequence of the equation of motion, the choreographic three bodies on 
the lemniscate satisfy the kinematical conservation laws independent of the potential 
energies, i.e., the center of mass (|^) and the angular momentum (^. Furthermore, they 
conserve a lot of kinematical and geometrical quantities, the kinetic energy (|TTp, the 



moment of inertia (|g), the product of mutual distance ([TsD, the sum of square of mutual 
distance ([l2|) and the sum of square of the curvature (p!0|). 

We have also pointed out that the tangent lines from the positions of the 
choreographic three bodies on the lemniscate meet at a point on the rectangular 
hyperbola, then have given the two geometrical methods to determine the choreographic 
three points on the lemniscate. We have given (i) how to determine the choreographic 
three bodies from a point on the rectangular hyperbola, and (ii) how to determine the 
choreographic three points from the given first point on the lemniscate. 

As we have already mentioned, we believed for some time the choreographic three 
bodies on the lemniscate was the Simo's figure eight motion because of the numerical 
similarity and a lot of conservation quantities of them. We hope that we can discuss in 
the future the orbit of Moore, Chenciner, Montgomery and Simo's figure eight motion 
in our context. 
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